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Through a detailed numerical investigation in three spatial dimensions, we demonstrate that 
long-lived time-dependent field configurations emerge dynamically during symmetry breaking in an 
expanding de Sitter spacetime. We investigate two situations: a single scalar field with a double- 
well potential and an SU(2) non-Abelian Higgs model. For the single scalar, we show that large- 
amplitude oscillon configurations emerge spontaneously and persist to contribute about 1.2% of the 
energy density of the universe. We also show that for a range of parameters, oscillon lifetimes are 
enhanced by the expansion and that this effect is a result of parametric resonance. For the SU(2) 
case, we see about 4% of the final energy density in oscillons. 



I. INTRODUCTION 



Spontaneous symmetry breaking plays a key role in our current understanding of particle physics and is expected 
to have been a major factor in determining the physical properties of the early universe [T]. In cosmology, two aspects 
of symmetry breaking are of great interest: it typically happens far from thermal equilibrium and it is inherently 
nonlinear. In the context of the electroweak phase transition, for example, an initially thermalized state is tossed out 
of equilibrium as the Higgs evolves to acquire a nonzero expectation value. In inflation, a nonthermal state thermalizes 
to reheat the universe with an explosive energy transfer from the inflaton to other field modes. It is thus of great 
interest to study the dynamics of symmetry breaking in an expanding background numerically in order to isolate key 
features that may escape analytical techniques. 

Here, we report results on 3d simulations for two situations: a single, self-interacting scalar field with a double-well 
potential, and an SU(2) non-Abelian Higgs model. In Refs. [3] results have been obtained for the case of a single 
scalar in Id. It was shown that long-lived, time-dependent field configurations known as oscillons j4j-[6] emerged 
spontaneously and contributed an amazing 50% of the total energy density. These initial results triggered the present 
study in the context of more realistic models. There are two broad classes of scalar field oscillons that have been 
studied in the literature, small and large-amplitude. Small-amplitude oscillons do not probe the highly nonlinear 
domain of the potential, and typically have large spatial widths [HJEIIE]. Their small amplitude makes it possible 
to study them using linearization techniques. Large-amplitude oscillons are harder to investigate analytically [5Hllj. 
Simulations of scalar models in static 2d and 3d backgrounds [12] and expanding Id backgrounds [5J [3] indicate that 
mostly large-amplitude oscillons are excited during symmetry breaking. As we show next, this is also the case for an 
expanding 3d spacetime. The situation is different for SU(2) models, as we explain below. 

This paper is organized as follows: in the next section, we introduce the scalar field model in an expanding universe 
and discuss its lattice implementation. We report our results for a double well potential, showing that oscillons 
contribute about 1.2% of the energy density. In section III we show that, contrary to naive expectation, for certain 
values of the expansion rate oscillons may have their lifetimes enhanced. We explain this result analytically by making 
use of parametric resonance. In section IV we introduce the SU(2) non-Abelian Higgs model and discuss its lattice 
implementation in an expanding universe. In section V we discuss the results for this model. In particular, we show 
that, as in the case of a real scalar field, oscillons contribute a nontrivial percentage of the total energy density. 
Furthermore, our results indicate that the cosmological expansion seems to favor the formation of oscillons for a wider 
range of parameters as compared to the static case, where oscillons were found only in a 2:1 mass ratio for the Higgs 
and gauge boson. In section VI, we briefly discuss possible application of oscillons in cosmology, which we hope to 
explore in forthcoming work, and conclude with a summary of our results. 
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II. SCALAR FIELD MODEL 



We consider a scalar field <&(x, i) propagating in (3 + l)-dimcnsional de Sitter spacetime with Hubble constant 
H = a/a and a double-well potential U(<&) — (A/4)[$ 2 — fi 2 /\] 2 . Using h = c = ks = 1 and defining dimensionless 
variables 4> — ^(m/V'-V) -1 and x v — [ix v (y — 0, 1, 2, 3), the equation of motion satisfied by (f> is 

a{t) a{t) 

where overdot and V denote derivatives with respect to dimensionless time x and space x l . The expansion rate 
becomes H = [iH, where H = d\n(a) / dx° . 

Our initial conditions simulate quasi-thermal states of the free massive scalar field. The parameters that control 
the distribution of the lattice modes are the temperature T and the mass of the field m = [i\/2. The simulation 
space consists of a cube with comoving size L and volume V = L 3 discrctizcd on a regular lattice with spacing Ax 1 — 
Ar (i = 1,2,3). We apply periodic boundary conditions and label the free field's normal modes by k = (27rnj/L), 
where n = (n x , n y ,n z ) and the rii are integers n$ = —A/2 + 1 . . . N/2. Here N — L/ Ar is the number of lattice points 
per side. Each free mode is described by a harmonic oscillator with frequency cj| = (2sin(/cAr/2)/Ar) 2 
k = |k|. The initial conditions for the field 4> are then given by 

*(r,t = 0) = ^E^K- ik - r + ^-' ik 1^ 

k 



r,t = = 



k 



where ak is a random complex variable with phase distributed uniformly on [0, 2ir) and magnitude drawn from a 
Gaussian distribution such that (|a/c| 2 ) = [coth(hu!k/2T) — l]/2. This is the amplitude distribution for a quantum 
harmonic oscillator [13] with the zero-point motion subtracted. On average, modes with hujk < T get assigned energy 
T, in agreement with equipartition, while the energy per mode goes rapidly to zero for fuvk > T. We thus need a 
lattice fine enough to resolve the high k modes that are excited at high temperatures. Using a value of Ar$ that is at 
least 10 times smaller than the wavelength of the mode satisfying huj^ ~ T is enough to provide a good continuum 
limit. 

We discretize the equation of motion using second-order space derivatives with lattice spacing Ar in all directions. 
We then step forward in time using a fourth-order Runge-Kutta method. By the Courant condition, we need to keep 
At < a(t)Ar at all times. We impose a maximum physical lattice spacing Ar max that is fine enough to resolve field 
configurations at physical sizes that we expect for oscillons. When a(t)Ar > Ar max , we refine the lattice by bringing 
the lattice spacing back to Ar max /2 and inserting points by polynomial interpolation. We pick Ar max and At small 
enough so that any further reduction does not significantly affect the final configuration of a run. All our simulations 
maintain energy conservation to a part in 10 3 or better. 

We evolve the field 0(x, t) in a box with 256 3 lattice points and Aro = 0.05/i^ 1 . We keep At = O.Ol/i -1 constant 
throughout the simulation. As the lattice spacing increases to Ar max = 0.5/^ _1 , we insert points in the lattice, bringing 
the spacing down to Ar max /2. We vary the values of the expansion factor H and the initial temperature T and evolve 
the field until it cools down to T/a(t) = 0.3/1. After the universe has expanded and cooled, we observed persistent 
localized structures as peaks in the energy density. We show a typical sequence of snapshots in Fig. [T] 

By isolating these peaks individually, we find that they all share the typical signatures of spherically-symmetric 
oscillon configurations: their centers oscillate with the typical oscillon frequency, as shown in the inset of Fig. [2] 
and their energies coincide with the plateau energies found in detailed oscillon studies [SJ |H1 HH] ■ The energy of a 
configuration is calculated by integrating the field's total energy around its peak using a radius r = lO/i -1 . We 
consistently found i? sc — 45/i/A. We then measured the fraction of energy in oscillons (Fig. [2]) and the number of 
oscillons nucleated as a function of temperature, which scales simply as N osc oc V — L 3 oc T 3 . 

Although we quote results for H = O.Ol/i, we have performed simulations for a slower expansion rate of H = 0.005/1 
obtaining similar qualitative behavior: for a wide range of initial temperatures, Posc/ptot = ^osc ~ 1.2%. Smaller 
values of H require impractical computation time, but we don't expect any qualitative changes. We note that since 
the simulations end with fairly large values of Ar = 0.5/Lt -1 , our results are lower bounds on p osc . Of course, scalar 
field oscillons are not stable in 3d and will decay after r osc ~ lO 4 /!^ 1 . Nevertheless, during their lifetime, they may be 
responsible for several important effects, as we discuss in section VI. We also note that oscillons are prevented from 
forming if the horizon size 1/H is of the order of the oscillon size, i? sc ~ 4/z _1 . For H > 0.1 fi, which fortunately is 
not very realistic, large-amplitude fluctuations are flattened out before the stabilizing effect of nonlinearities can kick 



FIG. 1: Sequence of time snapshots of the energy density. Time increases from left to right and top to bottom at times 
tfj, = 0,50,100,150,200,250, H — 0.01/i and T = 6.0/i. All snapshots show the energy density u with an isosurface at 
u = 0.2/A 
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FIG. 2: Fraction of energy in oscillons as a function of temperature in units of fi. Here Aro = 0.05/Lt - , Ar max = 0.5/x~ , and 
H — 0.01/i. Error bars denote ensemble averages over 10 runs. The inset shows the near-harmonic oscillations of the oscillon 
core. A simulation of a typical run can be viewed at http:/ /www.youtube.com/watch?v=_0co05XkNMY (web link on electronic 
version) . 



in. In other words, for oscillons to be cosmologically viable, we must have R osc /Xh — R OS cH <C 1, where Xjj = H 
is the horizon length. This condition is easily satisfied for physics below the Planck scale. 



III. LIFETIME ENHANCEMENT 



Having established that oscillons emerge dynamically in an expanding background, we need to examine how the 
expansion affects their lifetime. For numerical efficiency, we exploit the spherical symmetry of the final oscillon 
configuration and reduce our system to an effectively Id problem by letting V 2 — > d 2 (j>/dr 2 + {2/r)dcj>/dr in Eq. m 
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We find oscillons by setting the initial field configuration to be Gaussian, 0(r, 0) = 2exp(— r 2 / Rq) — 1, with boundary 
conditions cf>(r — > oo,<) = —1, <fi'(0,t) — 0, and cb(r, 0) = [5J In the absence of expansion, Gaussians with 
2.4 < Ron < 4.5 settle into long-lived oscillon configurations. 

We follow the same procedure as in 3d so that as soon as the lattice spacing becomes Ar max = O.f/i -1 , L max > 2/H, 
we insert points via polynomial interpolation, and bring the lattice spacing back to Ar = 0.05/i _1 . We then truncate 
the box to L > 1/H, which can't affect the oscillon at r = 0. We always use a box of initial size Lq = 1/H + 
in natural units, and we have verified that any run with Lq > 1/H gives identical results. In Fig. we show the 
effects of expansion for a sample of initial configurations. There is a clear symmetry about Rq — 2.86/x , the longest- 
lived oscillon in the absence of expansion: radii to both sides of R = 2.86/i" 1 experience an increase in lifetime 
for a range of H, with the increase being more pronounced for shorter lifetimes. The longest-lived oscillon, in turn, 
doesn't experience any noticeable enhancement. The inset of Fig. [3] shows the maximum fractional increase in lifetime 
(T max — To) /to as a function of initial radius Rq. The lifetime enhancement follows an approximate scaling law around 
R Q = 2.86/i-\ r oscf i ~ \R of i - 2.86| 05 . 




FIG. 3: Lifetime, in units of jj, , for oscillons formed from Gaussians with 2.4 < R/j, < 4.4, as a function of expansion rate H 
in units of fi. The inset shows the maximal fractional increase in oscillon lifetime for different radii at i/ max . 

To understand the origin of the lifetime enhancement caused by the expansion, we decompose the field as </>(x, t) = 
4>av(t) + t), where <f> av is the volume averaged field. Linearizing Eq. [I] with respect to S(f>(ic,t) and taking the 
Fourier transform, we obtain (for k > 0) 

5$(k,t) + 3H8i(k,t) + (J^ + V" (<t> av (t)))ty(M) = 0. (3) 

Once the Gaussian has settled into the oscillon stage, V"((f> av (t)) can be approximated by V"((p a v{t)) = $o cos(W) + C, 
where $o and C vary very slowly during an oscillon's lifetime and depend on the value of the initial radius Rq and 
the Hubble constant H. Here to < m is the oscillon's frequency of oscillation. Introducing new variables u>t = 2z — n, 
and 5(f) = exp(—3Hz/uj)x, Eq. [3] becomes 

X " + [A k -2qcos2z} X = 0, (4) 

where A k = ± [4k 2 /a 2 + AC - 9H 2 ] , q = 2$ /^ 2 , and prime denotes differentiation with respect to the new variable 
z. Eq. [i] is the Mathieu equation, which is known to exhibit parametric resonance when Af. ~ I 2 , I = 1,2,.. . [14] - 
Thus, particular combinations of values of C, H and u can lead to exponential amplification in the oscillations of X) 
and consequently 8<j)(k, t), at certain modes k/a. Because C is a positive number, the A k — 1 resonance window occurs 
for real values of k/a only for very large values of H, which destabilize the oscillon before it can lead to resonance. 
The A). ~ 9 and higher windows lead to resonances that are too weak to overcome the damping due to the expansion. 
The A k — 4 window, however, can lead to parametric amplification of the dominant oscillon wavevectors for the values 
of H that generate the observed lifetime enhancement depicted in Fig. [3] 

The Ak — 4 resonance window leads to exponential amplification to the oscillations in x oc e^ 2 , where £ ~ V5q 2 /4:8 
[15] . For small values of if, the amplification overcomes the damping due to expansion. For higher values of H, the 
damping overcomes the amplification, and the oscillon decays. An example is shown in Fig. [4j 
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FIG. 4: Oscillations in S(f>(k,t) for an oscillon with Ro = 2.5/i -1 and expansion rate H = 0.006/t, at times 3200 < tjj, < 3700. 
Here we plot a mode which starts at k/a ~ 1.4/1 and gets redshifted to k/a ~ 0.1/i. As k/a gets smaller, the mode enters the 
resonance window, and its Fourier component gets amplified and then redshifted away. In this case, ui ~ 1.4/i (< m = \/2/t) 
and C = 1.8/t 2 , with Ak = 4 for fc/a ~ 0.4/i, which is the dominant wave vector of this oscillon (R ~ 2.9/i _1 ). 



IV. SU(2) MODEL: IMPLEMENTATION 



Because the SU(2) model is considerably more expensive numerically to simulate, the range of experiments we can 
carry out is limited. However, these experiments show very similar behavior to the scalar model. We begin from the 
Lagrangian density in the absence of expansion, 



C 



1 



2\2 



F^ ■ F"" + (ZJ /1 $)t(X)f $) _ A(|$| 2 - 



(5) 



where the boldface vector notation refers to isovectors. Here the Higgs field $ is an SU(2) doublet, and the SU(2) 
field strength and covariant derivatives are 



F„ v = d^W v - d v W^ - gW^ x W v , 

D^F^ = WF^ - gW» x F M „ , 
where r represents the weak isospin Pauli matrices. We obtain the equations of motion 

D^F> XV = J V , D^D^ = 2A(/i 2 - |$| 2 )$, 



(6) 
(7) 
(8) 

(9) 



where the gauge current is J v = glm (Z?„$)^t<1> and we work in the gauge Wo = 0. With this choice, the covariant 
time derivatives become ordinary derivatives and we can apply a Hamiltonian formalism. The Wj fields have mass 
m w = ff/Vv^, and the Higgs field has mass = 2[i\f\. 

To include the effects of expansion, we again work in comoving coordinates with a scale factor a(t). We now have 
the action 



S 



d 3 ra{t) 3 



3 B 3> 



where 



E j =W j and B j 



j'J"=x,y,z 

Here dot indicates time derivative and Latin indices run over space dimensions 



(10) 



(11) 



G 



For numerical computation we put the theory on a lattice, following the conventions and techniques used in 116) . 
The field variables are the values of the $ p field at the lattice sites p and the spacelike Wilson lines 



U P (t) = e l 9Wj(t)-Ta(t)Ax/2 



(12) 



emanating from lattice site p in the spacelike direction j. Since the lattice equations are second order, we will find 
each field at the next time slice based on the previous two. We let t be the time for the current set of lattice points 
and spacelikc links and define t + = t + At and t_ = t — At to be the subsequent and previous times respectively. 
We also take t + / 2 = t + At/2 and t_/ 2 =t — At/2 to be the times in between, which will be the times at which we 
evaluate the timelike links. 

We define the Wilson line for the link emanating from lattice site p in the negative j th direction to be the adjoint 
of the corresponding Wilson line emanating in the positive direction from the neighboring site, U^_j(t) = Uj~ J (ty, 
where the notation p± j indicates the adjacent lattice site to p, displaced from p in direction ±j. At the edges of the 
lattice we use periodic boundary conditions. We define the elements of the field strength tensor, which are centered 
on the timelike and spacelike plaquettes of the lattice, 



E p 3 {t +n 



iga(t 



+/2 



}AxAt 



\ogU P (t+)U P (rf and r-B p At) 



g(a(t)Ax) 2 E ' 

j',j"=x,y,z 



U'rVutf ( 13 ) 



where U^ Ujl) (t) = U P {t)U P + 3 (t)U p + 3+3 \t)U p ^/ {t) and we have defined the logarithm of a 2 x 2 matrix in the form 



of Eq. 12 



as 



iga(t)Ax 



\ogU p {t) = ^f^W p {t) 



(14) 



We note that \ogXY ^ logX + logF when the matrices do not commute. The logarithms and exponentials needed 
to convert between the group and the algebra can be computed efficiently using 



in ■ t sin a = 



cos 9 + in z sin v in x sin v + n v sn 
n x sin Q — fiy sin 6 cos 9 — in z sin ( 



(15) 



where h is a unit vector and the link matrices have fiQ = W p (t)ga(t) Ax/2. For efficiency we replace sin# — > 9 and 
cos 9 — > yl — 9 2 when computing both the logarithm and the corresponding exponential. This discretization then is 
equivalent (without expansion) to what is used in other numerical studies of electroweak symmetry breaking dynamics 

We find the equation of motion for the Higgs field 



3HAt 



3HAt\ 



2$ p (t) - I 1 — J $ p (t_) + At 2 $ p (t) 



(16) 



where H = is the Hubble constant and 



$ p (t) 

For the gauge fields, we have 



E 



U p {t)<S> p+j {t) - 
a{t) 2 Ax 2 



2A - |$ p (t)| 2 ) $ p (<) . 



(17) 



U p (t + ) = (expjlog^Ci)^)^- 1 ) 



v / to8^ ) (t)+logP5 (3 -,- J -> ) (*) \ ia(t)Ax 
^ 1 a(t) 2 Ax 2 I + 9J i W ' 



where the gauge current is 



J p {t)=glm 



At 2 ) U p {t) 



^ p {t)^TU P {t)^+At) 

a{t)Ax 



(18) 



(19) 
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and the logarithm in Eq. [14] is used to compute the exponents in Eq. 18 

Assuming it is obeyed by the initial conditions, time evolution preserves the Gauss's Law constraint, 



E 



iff(t +/2 ) + (t +/2 
a{t + / 2 )Ax 



(20) 



where the charge density is given by 



Ja(t+/2) = 5 Im 



$ p (i+) - $ p (t) 
At 



T^> P (t) . 



(21) 



Energy is not conserved because in the expanding background we have dll — —pdV, where p is the pressure. In the 
lattice model we then have 



dU 
lit 



-(a(t)Ax) 3 Hj2 



\ E & 



-■El 



BIB-) 



3|$l 



E 



a(t)Ax 



3A(|$ P | 2 -M 2 ) 2 



(22) 



We have verified that both the Gauss's Law and energy constraints are well obeyed throughout our simulation. 

We set initial conditions for the first two time-slices, which we denote as to and ti, by occupying the modes of all 
the components of the Higgs and gauge fields at temperature T, as in the case of a single scalar field. Then we modify 
these initial conditions to make them obey Gauss's Law, by the following steps: 

• First, since we have periodic boundary conditions, the total charge should be zero. To enforce this constraint, 
we shift the <& p field on both of the first two time slices to and t\ by the same constant, 



$p $p 



(Jo • t) 4 , 



(23) 



where Jo and $ are the average values of Jo an d 4 over the lattice at time tin, respectively. 

• Next, we fix the longitudinal component of the gauge fields, as described in Ref. [18] . We take a discrete Fourier 
transform of the initial charge Jq (£1/2) and gauge field Jj'(*i/2) to obtain Jq (^1/2) an d Wf(t 1 / 2 ) for the initial 
time step, where k labels the Fourier transformed lattice. We then modify the initial time derivative of Wj (by 
changing its value on one of the first time slices but not the other) by sending 



w](t l/2 ) -+ w](t 1/2 ) - ( **'WV(ti/ 2 ) + i4{h,2) 



n 3 



(24) 



and then inverting the discrete Fourier transform to obtain the modified gauge fields, which in turn give the 
modified Wilson loops. Note that we don't make any modification for k = 0, where this transformation breaks 
down; that case was already handled by the previous step. 

• Finally, while in an Abelian theory the subtraction of the longitudinal component would be sufficient to imple- 
ment Gauss's Law, for a nonabelian theory the nonlinear term in the field strength makes this agreement only 
approximate. As a result, we adjust the phase of $ p (ti), 



\$ p (h)\ 



W(t 1/2 )*P(t ) with W{t l/2 ) 



exp 



E 

j=x,y,z 



logC/f(i 1 )C/f(io) t + log^,(i 1 )^,(i )t 



5 2 (a(t 1/2 )Ax) 2 |<I>f(i 1 )||^(io)|/2 



leaving $q(*o) unchanged, in order to assure that Gauss's Law is satisfied. 



(25) 



V. SU(2) MODEL: RESULTS 



We begin with a universe of size La(t = 0) = 4//i, temperature T = 4/i, lattice spacing a(t = 0)Ax = 1/(224//), 
and use a time step At — 1/(448//). We allow the universe to expand at a constant rate, with Hubble constant 



H = /i(log2)/12 ss O.O6/1, and expand the universe by a factor of 224, so that the final lattice spacing is a(tn na i)Ax = 
1/ fx. We measure the fraction of energy in oscillons by including those points whose energy density is four times the 
average energy density. Under ordinary thermal expansion this fraction would stay constant, and during the initial 
stages of the expansion it is identically zero. At the end of the expansion, approximately 4% of the energy is found in 
oscillons by this measure. The energy density is shown in Fig. [5] and the evolution of the energy over time is shown in 
Fig. [6] Compared to the scalar model, oscillons in the SU(2) model have smaller amplitude and larger spatial extent, 
requiring a larger simulation volume. The simulation is thus considerably more expensive numerically, especially given 
the cost of evolving a total of thirteen real degrees of freedom per lattice site instead of one. 




FIG. 5: Energy density in units of fj, 4 at the end of the expansion for g — y2 with A = 1.0 (left panel) and A = 1.1 (right 
panel). To implement the Standard Model coupling of gsM = 0.624 with the same mass ratio, this energy would be scaled up 
by a factor of (g/gsyi) 2 - 
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FIG. 6: Fraction of energy in oscillons (solid line) and total energy (dashed line) as functions of time. Energy is given in units 
of fi and time in units of l//i, with g = \/2 and A = 1.0 (left panel) and A = 1.1 (right panel). To implement the Standard 
Model coupling of gsM = 0.624 with the same mass ratio, this energy would be scaled up by a factor of (<?/</sm) 2 - 



We note that similar results are seen both for g = \/2, A = 1, in which the Higgs and gauge fields are in the 
2 : 1 mass ratio found in [TBI 123], and for g = \/2, A = 1.1, where the masses are not in this ratio. Because of the 
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high numerical costs associated with the expanding background simulation, we are not able to track the stability of 
oscillons formed in this way over long time scales, but these results suggest that the expansion may broaden the range 
of parameters for which oscillons are stable. Also, by allowing oscillons to form from a thermal background rather 
than a fixed ansatz, this simulation is capable of scanning a wider range of configuration space (and our results clearly 
show that oscillon configurations represent attractors in this space). Work is currently underway to investigate these 
questions in greater detail. 

VI. POSSIBLE IMPACT ON COSMOLOGY AND SUMMARY 

Our results indicate that oscillon-like configurations emerge dynamically during spontaneous symmetry breaking in 
expanding cosmological backgrounds. Furthermore, they contribute a significant fraction of the total energy density. 
They are thus poised to play an essential role in the dynamics of the early universe, be it during post inflationary 
reheating or during symmetry-breaking phase transitions. 

In order to briefly address the impact oscillons may have on cosmology, it is best to consider different energy scales 
separately. For the sake of illustration, we focus on the GUT and electroweak scales, which differ by roughly 13 orders 
of magnitude. Also, it is important to differentiate between real scalars and Abelian and non-Abelian Higgs models. 
Thus, before we start, it may be useful to summarize what is known of oscillon lifetimes in these models. 

As we mentioned before, for real scalar fields in 3d, the oscillon lifetime-with or without the enhancement from the 
expansion reported here-is typically t osc ~ 10 3 ~ 4 /i -1 [SJ[TU]. For models with gauge fields, the evidence at hand points 
to very large lifetimes. Studies for Abelian-Higgs models in 2d have not seen oscillons decaying, and report lifetimes 
in excess of lO 5 /!^ 1 [21] . Studies of Abelian-Higgs models in 3d obtained similar results: oscillons have been observed 
to persist for times t > 7 x lO 5 /!^ 1 without decaying [23. For non-Abelian Higgs models, the situation is similar: 
the data at hand indicates that once formed, oscillons live for extremely long times. Those in the gauged-SU(2) 
Higgs model have not been observed to decay after t > 5 x 10 5 // -1 [TO1 [231 . Thus, although a more detailed study 
of Abelian and non-Abelian Higgs oscillons and their stability is clearly warranted, results so far indicate that they 
may be extremely long-lived, even perturbatively stable. Of course, the key question is whether their lifetime can be 
longer than the cosmological time scale at their formation. If that's the case, they behave as stable, localized defects. 

At the GUT scale, it is clear that the oscillon lifetime is at least of order of the cosmological time scale, H^ 1 ~ 
(-Mpianck/A 4 ) ~ 10 3-4 . Thus, for all practical purposes at GUT scales oscillons behave as stable localized defects. As 
has been shown elsewhere, in the context of first-order phase transitions, long-lived bubble-like configurations such as 
oscillons can either become a critical bubble or coalesce to become one. In both cases, the decay of the false vacuum 
is greatly accelerated, changing from exponentially-suppressed to power-law [24| . It has been suggested that oscillons 
may accelerate the decay of the false vacuum during inflation, potentially solving the bubble coalescence problem of 
old inflation. This has been recently illustrated within the context of a modified hybrid inflation model |25j . 

Oscillons may also have a key impact during post-inflationary reheating. As coherent field configurations, they 
naturally delay the approach to equilibrium, acting as bottlenecks for equipartition [26j . As such, they may influence 
(decrease) the reheating temperature, a possibility we are currently investigating. An interesting open question is 
how these nonequilibrium results apply in the context of gauge models. 

Moving on to the electroweak scale, since H ew ~ 10~ 16 , we are on a realm which is very distant from our numerical 
range of H ~ 10~ 2 . Still, we suggest that there are at least two ways in which oscillons may play a role at these 
relatively low energy scales. Both depend on their lifetime. If non-Abelian oscillons live for t ~ lO 16 /! -1 , that is, if 
they are perturbatively stable, they will remain active at cosmologically-relevant time scales. As at the GUT scale 
discussed above, they may speed up vacuum decay in the context of a first-order transition (which is ruled out in the 
Standard Model but not in all of its extensions) or they may delay thermalization. If they persist for even longer, 
they may even be relevant to dark matter or baryogenesis. 

On the hand, if they live for shorter times 10 4 < t[i < 10 16 , their presence may still affect the dynamics of symmetry 
breaking. As is well-known, most phase transitions are initiated due to the presence of inhomogeneities or "seeds" [13] . 
If oscillons are present in sufficient quantities, they will modify the effective potential nonperturbatively, affecting the 
dynamics of the transition [37]. Although much work remains to be done to investigate such nonperturbative effects 
in more detail, these mechanisms suggest that even relatively short-lived oscillons will have important effects during 
cosmological symmetry breaking. 
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Note added in proof: A recent work by Mustafa A. Amin offers further support to our hypothesis that oscillons will 
have important effects in an expanding universe |28j . We thank the author for sending his manuscript to us. 
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